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Abstract. Integrable discretizations are introduced for the rational and hy- 
perbolic spin Ruijsenaars-Schneider models. These discrete dynamical systems are 
demonstrated to belong to the same integrable hierarchies as their continuous-time 
counterparts. Explicit solutions are obtained for arbitrary flows of the hierarchies, 
including the discrete time ones. 



1 Introduction 



The famous Calogero-Moser many-particle model possesses the most rami- 
fied tree of all possible kinds of generalizations. Even at the classical (non- 
quantum) level, one has three types of the usual Calogero-Moser (CM) sys- 
tems: rational, hyperbolic/trigonometric, and elliptic [HI, @. They may be 
viewed as being related to the Ai^_i root system, and there exist immedi- 
ate generalizations to other root systems [Q. A relativistic generalization 
of the CM systems constitute the Ruijsenaars-Schneider (RS) models [Q], 
[^], appearing also in three variants (rational, hyperbolic/trigonometric, and 
elliptic). Some time ago spin generalizations of the CM models were intro- 
duced 0, 0, recently they attracted more attention |^, 0, 0. Finally, a 
spin generalization of the RS model was introduced and studied very recently 
[p!0[] ; it found its application also to the solitons of affine Toda theories 



(We have not mentioned some other directions of generalizing the CM model, 
such as imposing external potentials etc.) 

Still another direction of the development of CM type of models is their 
time- discretization, which was started by finding the discrete time CM system 
in [|12| and followed by the discretization of the RS system in |jl3|. (Actually, 
the discretizations of the CM model and some its generalizations appeared a 



decade earlier in the image of Backlund transformations [14]). The integrable 



discretizations found in these papers turned out to have a remarkable prop- 
erty: they belong to the same integrable hierarchies as their continuous time 
counterparts, i.e. they share with continuous time systems such attributes as 
Hamiltonian structure. Lax matrices, r-matrices, action-angle variables and 
so on. This property was clearly identified and put in a basis of a general 
procedure of discretization of integrable systems in [jT5[, [|T^, for the further 



development see [[Tj, [|T8|, pJ [. 

In the present paper we continue this line of research by applying this 
general procedure to the spin RS model. Despite the fact that a Hamiltonian 
structure for this recently introduced system is still unknown, it is possible to 
attach a hierarchy of fiows to this model and to solve explicitly an initial value 
problem for a general fiow of the hierarchy. This is done (for the rational 
and hyperbolic systems) in Sect. 3. The discrete time systems are introduced 
and studied in Sect. 4. In Sect. 5 we show what do our results give for the 
usual RS model (without spins): even in this case we obtain new results, or 
better, give new insights. Several open problems are outlined in Sect. 6. 
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2 Spin RS model 



We prefer to give here the results for the hyperbohc and rational spin RS 
models in a form suitable for our present purposes, since it is difficult to 



extract the corresponding results from the general setting of [jTO . 

The model is described in terms of scalar variables Xk, 1 < k < N, N 
vectors \ak) G R"', 1 < A; < A^, and covectors {bk\ G R'^, 1 < A; < A^. We 
shall represent vectors \ak) as columns of the d x N matrix A, and covectors 



->k\ 



as rows of the N x d matrix B. 



The spin RS model proper as given in []I0| is described by the equations 
of motion: 

ik = {bk\ak), (2.1) 
\dk) = -XkWk) -Yj\^j)i^j\^k) V{xj - Xk); (2.2) 

{bk\ = Xk{bk\ + '^{bk\aj){bj\V{xk-Xj). (2.3) 
Here for the rational model 

11 T 

V{x) = = , ^ (2.4) 

X X + 7 x{x + 7) 

and for the hyperbolic one 

sinhfo') 

V(x) = coth(a;) - coth(x + 7) = — ri—^ — tt 7, (2.5) 

^ ^ ^ ^ ^ ^ smh(x)smh(x + 7) ^ ^ 

7 being parameter of the model. Further, Xk = \k{x,a,b) in (|2.2|), ( |2.3| ) 
are arbitrary functions of the dynamical variables. The simplest possible 
choice is, of course, Afc = 0. Solutions of the system with arbitrary A^'s may 
be obtained from the solutions of the system with = upon change of 
variables \ak) ^ Ak\ak), {bk\ ^ A^^{bk\ with = exp Xk{t')dt'^. 
These models admit Lax representations 

L = [M, L] (2.6) 

with the Lax matrix, which we prefer to choose in the form 

^ 7 

L = L{x,a,b) = — —{bj\ak)Ejk, (2.7) 

j,k=l -^i ^k + 1 
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or 

L = L{x,a,b) = f: -^^^^^^^—^AbMk)E,,, (2.8) 

in the rational and hyperbolic cases, respectively. The matrix M = M{x, a, b) 
has the entries 

Mjk = Vixj - Xk){bj\ak), J ^ k, (2.9) 

Mkk = \k. (2.10) 

In fact, provided the ansatz ( |2.7|) or (|2^ ) is given, the Lax equation (|3.6| ) 
may be derived as a consequence of the following two equations, which are 
in a sense more fundamental: 

A = -AM, B = MB. (2.11) 

So far, to our knowledge, the Hamiltonian formulation of this model has 
not been derived. This makes it impossible to give a rigorous statement 
about complete integrability of the spin RS model. However, an exphcit 
theta-function solution for this model was given in the [1^, and we shall 



now demonstrate that in the rational and hyperbolic cases a simpler solution 
can be given, in terms of eigenvalues of certain matrices constructed explicitly 
from the initial conditions, just as for usual CM and RS models ||l|], 0, and 
for the spin CM model ||^. 



3 Spin RS hierarchy 

We want to show now that the rational and the hyperbolic spin RS mod- 
els described above are in fact only the simplest representatives of the two 
corresponding hierarchies, and to derive Lax representations for each flow 
of these hierarchies. The members of the hierarchies are encoded by an 
arbitrary conj ugation-cof anant function f{L) on gl{N), the spin RS mod- 
els proper corresponding to the case f{L) = L. (Note that in the hypo- 
thetical Hamiltonian formulation the members of the hierarchies would be 
encoded by conjugation-mf anani Hamiltonian functions <f{L), such that 
f{L) = LVipi^L) = V(p{L)L, cf. in the non-spin case.) 
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For a conjugation-covariant function f{L) we introduce the system of 
differential equations: 

Xk = f{L)kk, (3.1) 
W) = -Afc|afc) - \aj)Mjk, (3.2) 

{bk\ = Xk{bk\ + Y.Mkj{bj\; (3.3) 
here in the rational case L is defined by (p.7|), and 

M,k = ^^. J^k; (3.4) 
in the hyperbolic case L is defined by ( |2.8|) , and 

M^k= . J {L)jk, J^k- (3.5) 

sinh(a;j — Xk) 

the scalars A^, 1 < A; < are, as before, arbitrary functions on x, a, h. (It is 
easy to see that we recover the system (|2.1| )- ([2l^ ) for f{L) = L). 

As the spin RS model proper, this general system also admits Lax repre- 
sentation 

L=[M,L], (3.6) 
which, given the ansatz (|2.7|) or (p^.S]) for the Lax matrix, is a consequence of 



A = -AM, B = MB. (3.7) 

Here the off-diagonal entries of the matrix M are given by ( p.4| ) or ( |3.5| ), 
respectively, while the diagonal ones - by 

Mkk = Xk. (3.8) 

To solve the system explicitly for an arbitrary f{L), we need to introduce 
an auxiliary diagonal matrix X by the following formulas: in the rational 
case 

X = diag(a;i, . . . ,X7v), (3.9) 

and in the hyperbolic case 

X = diag(e2^S...,e2^'^). (3.10) 
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It is important to note that the ansatz ( p.7| ) for the rational Lax matrix 
is equivalent to 

XL- LX + -fL = jBA, (3.11) 

while the ansatz (|2.8| ) for the hyperbolic Lax matrix is equivalent to 

e^XLX-^ - e-^L = 2smh{-f)BA. (3.12) 

A key observation necessary for an explicit integration of equations of 
motion, is the following: the evolution equation for Xk ( |3.1| ) together with the 
expressions ( |3.4|) , (|3.5| ) may be cast in a single matrix differential equation for 
X, whose form depends on the case under scrutiny. Namely, for the rational 
one: 

X=[M,X] + fiL), (3.13) 

and for the hyperbolic one: 

X =[M,X] + 2Xf{L). (3.14) 

Now we proceed as in |l[], 0, p: define the matrix V = V{t) as the 
solution of the matrix differential equation 

VV-^ = M (3.15) 

with the initial condition V"(0) = /. Then in the rational case 

X = X{t) = V{X, + tf{L^))V~\ (3.16) 

in the hyperbolic case 

X = X{t) = VXQe^'^'-^°'^V-\ (3.17) 

and in the both cases 

L = L{t) = VLoV-\ (3.18) 

A = A{t) = AoV-\ B = B{t) = VBo. (3.19) 

So, just as in the case of the usual CM and RS models and the spin 
CM model, in the rational case the positions are the eigenvalues of the 
matrix Xq + tf^Lo), and in the hyperbolic case exp(2xk) are the eigenvalues 
of the matrix Xq exp(2t/(Lo)). As in the spin CM model, the corresponding 
eigenvectors (i.e. the diagonalizing matrix V) define the evolution of spin 
degrees of freedom. The non-uniqueness of the diagonalizing matrix (which 
is in principle defined up to a left multiplication by a diagonal matrix) is 
reflected in the arbitrariness of the diagonal part of the matrix M = VV~'^, 
i.e. of the scalars (cf. (|378|)). 
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4 Integrable discretizations of the spin RS 
models 



We introduce now two families of maps 

(xfc, \ak), {bk\) ^ (xfc, iSfc), (bk\), 

depending on the small parameter h > and approximating, as h —>■ 0, the 
rational and hyperbolic spin RS models, respectively. Just as the spin RS 
models have N free parameters = Xk{x,a,b), will our maps have N free 
parameters /i^. 

The rational spin RS map is defined by: 

Xk - Xk = hfi^^(bk\ak), (4.1) 
- laO = -h { E - E '"^^''^'^^ ] , (4.2) 

yjyfc ^3 -"^k j=l Xk + 1 J 

{K\-,m=hh^^^^-j:B^], (4.3) 

^fc j=i -^fc ^ 1 j 

where the parameters [ik are supposed to satisfy the following asymptotic 
relations: 

/ifc(x, x, a, a, 6, 6; /i) = 1 + {x^ x, a, a, 6, 6) + 0{h^). (4.4) 

The simplest possible choice is, of course, /.t^ = 1. 

The standard implicit functions theorem implies that the equations ( [4.1|) - 
(^4.3|) for h small enough have a unique solution (x, a, h) which is 0(/i)-close to 
(x, a, 6). Hence a map (x, a, h) ^ (x, a, 6) is defined, being an approximation 
of the time h shift along the trajectories of the rational version of the system 
(|3)-(PD with 

Xk{x,a,b) = n'j^\x,x,a,a,b,b) - 7"^(&fc|afc) (4.5) 

(the second term in the right-hand side of this formula stems from the fact 
that in the second sums in (|4.2| ), (|4.3| ) the term corresponding to j = A; is 
present, as opposed to (|2.2|), (P^). The choice /ifc = 1 corresponds, therefore. 
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to the flow of tlie rational spin RS model with = —7 ^{hk\ak). In order 
to get an approximation for the flow with = one can take, for example, 
/ifc = 1 + /i7-^(6fe|afc). 

The hyperbolic spin RS map is defined by: 



g2{a;fc-Xfc) _ 



-h 



\aj){hj\ak){cot\i{xj - Xk) 



N 



\aj){hj\ak){cot\i{xj - + 7) - 1 



(4.6) 



(4.7) 



{bk\ - f^kihl 



h {bklo-i) {bj I (coth(xfc — Xj) — 1 

N ^ 

-J2(^k\aj) {bj\(coth.{xk - Xj +7) 



(4.8) 



Under the same assumtion ( [4.4|) as before, this system defines a map (x, a, b) h-> 
(x, a, b) approximating the time h shift along the trajectories of the hyper- 
bolic system (|2.1|) - (p.3|) with 

Xk{x,a,b) = fi'j^\x,x,a,a,b,b) - (coth(7) - l){bk\ak). (4.9) 

The choice Hk = 1 corresponds this time to the flow of the hyperbolic spin 
RS model with A^ = — (coth(7) — l){bk\ak). 

An outstanding property of the introduced discretizations for the spin 
RS models is their solvability. More precisely, they turn out to belong to 
the same hierarchies as their continuous time counterparts. That means that 
they admit discrete time Lax representations with the same Lax matrices as 
the continuous systems. 

Theorem 1. The maps (^ )-(^) and (^) - (^) have the following 
discrete Lax representations : 



LM =ML 
AM = A{I + haL), {I + haL)B = MB. 
Here for the rational case a = 7"^, L is as in (|2.7|), and 

h{bj\ak) _ 



M. 



jk 



Xj Xjz 

7 



(4.10) 
(4.11) 

(4.12) 



for the hyperbolic case a = coth(7) — 1, L is as in ( p.8| ), and 



M 



jk 



sinh(xj — Xk) 



{bj\ak); 



in the both cases the condition 



■M-kk — fJ-k 



(4.13) 



(4.14) 



is imposed. 

The proof of this Theorem follows by direct calculation. It should be 
remarked that the same phenomenon as in the continuous cases takes place, 
namely that the evolution equations ( [4.11|) for A, B are in a sense more 
fundamental then the Lax equation ([4.10| ) for L: given the ansatz ( |2.7| ) or 
(|2.8| ) for L and a corresponding ansatz ( |4.12| ) or ( [4.13| ) for A^, the Lax 
equation ( [4.10| ) is a consequence of ( ^.11| ) 



From the Lax equation ( [l.lOj ) follows that the equations ( 4.11 ) can be 
presented also as 

A = A{I + haL)M-^, B = M{I + haL)-^B. (4.15) 

As a consequence we get: AB = AB, or 

N ^ N 

j=i j=i 

This allows to rewrite the hyperbolic map ( [4.6|) -( ^^) in the following form: 



tanh(xfc — Xk) 
^k\ak) - \ak) 



hi^k {bk\ak), 

-h( \aj)(bj\ak)coth.{xj - Xk) 

N 

- \aj) {bj\ak) coth{xj - + 7) 
i=i 



(4.16) 



(4.17) 



{bk\ - i^k{h\ = hi Y{bk\aj){bj\coth.{xk - 



N 



Y{bk\aj) {bj\coth.{xk - Xj +7) 
i=i 



(4.18) 
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(where z/^ = /Xfc + h{bk\ak), so that, for example, the choice t'fc = 1 leads to a 
discretization of the hyperbolic spin RS model with = — coth(7)(6fc|afc))- 

It turns out to be possible to find continuous time flows from the spin RS 
hierarchies interpolating our maps (we say that a map is interpolated by a 
flow, if each orbit of the map is given by the values at t = nh, n G Z, of a 
certain orbit of the flow). The following statement gives these interpolating 
flows. This, further, provides us with explicit solutions of the introduced 
discrete time dynamical systems. 

Theorem 2. The map ([4.1|) - ([4.3|) is interpolated by the flow of the ra- 
tional spin RS hierarchy ( p.l| )-( pl3| ) with the function 



1 r\-l 



fiL)=LiI + h^-'L) 

Analogously, the map ( [4.6|) -( ^l8D is interpolated by the flow of the hyperbolic 
spin RS hierarchy ( ^.1D -( PT3D with the function 

/ j + /i(coth(7) + l)L \ 
■'^ ' 2/i l^/ + /i(coth(7) - l)Ly ■ 



Proof for the rational case. Equations of motion ( 4.101 ), ( [4.15 ), where 



a = 7 , have to be supplemented with an ansatz for A4, which may be also 
represented as an equation of motion for the auxiliary diagonal matrix X 
from (pD: _ 

XM - MX = KB A. 

Introduce now the matrix 

M = M{I + h-i-^ L)-^ . 

In terms of this matrix evolution equations may be presented as 

L = MLM-\ (4.19) 

A = AM-\ B = MB, (4.20) 
XM{I + h-f-'^L) - M{I + h-f-^L)X = HMBA. (4.21) 
But the ansatz (|2.7|) for L, being equivalent to the commutation relation 



(|3.11|) , implies: 

(/ + h-f-^L)X = X{I + /i7"^L) + hL- KB A. 
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This allows to rewrite ( |4.21| ) in the form 

XM - MX = hML{I + h-i^^L)-^ = hMf{L). (4.22) 

Defining now V as the solution of the matrix difference equation 

t/y-i = M (4.23) 

with an initial condition Vq = /, we see immediately that the solution of 
( gl9D , iWM) . O) is given by 

L = VLqV'\ A = AoV-\ B = VBo, (4.24) 

X = V{Xo + nhf{Lo))V-\ (4.25) 

which proves the Theorem in the rational case. 

Proof for the hyperbolic case. The equations of motion (f4.10 ), ( 4.15|) , 
where this time a = e^'''/sinh(7), have to be supplemented with an ansatz 
for Ai , which is this time equivalent to the following equation of motion for 
the auxiliary diagonal matrix X from ( |3.10|) : 

XM - MX = 2hBAX. 

In terms of the matrix 

M = M{I + haL)-^ 

the evolution equations for L, A, B take the form ( [4.19| ), ( |4.20|) , and for X 
- the form 

XM{I + haL) - M{I + haL)X = 2hMBAX. (4.26) 

The ansatz ( |2.8| ) for L is equivalent to the commutation relation ( |3.12| ), 
which, in turn, implies: 

(/ + haL)X = X{I + hae^^L) - 2hBAX 

(where we have taken into account that ae'''sinh(7) = 1). This allows to 
rewrite ( [4. 261 ) in the form 

XM = MX^-^^^^^ = MXe^^f^'^l (4.27) 
/ + haL ^ ' 

Defining now V by ( |4.23| ), we see immediately that the solution of (^4.19| ), 
(CT) , (iH is given by (1041 ) and 

X = VxA :[±i^^^y y-i ^ vx,e^^^^^'^^^V-\ (4.28) 
which proves the Theorem in the hyperbolic case. 
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5 Connection with the usual RS model 



In the case d = 1, i.e when \ak), {bk\ become scalars, the spin RS models 
yield back the usual spinless ones @ , @] . The usual RS hierarchy is now well 
studied. In particular, explicit solutions were found in 0, for the rational 
and hyperbolic cases, and in |jlO| in the elliptic one; the Hamiltonian structure 
of the usual RS hierarchy is known and admits an r-matrix interpretation 

n, n, m, n. 

To obtain the usual RS model from its spin counterpart in the case d = 1, 
one sets a^bk = Ck in (|3.1| )- (|3.3| ), and it follows irrespective of the values of 
the parameters Xk'- 

Ck = CkY^ Cj [V{xk - Xj) - V{xj - Xk) . 

These equations describe the usual RS model, the simplest flow of the usual 
RS hierarchy. 

The same reduction is admissible for our discrete systems. Namely, ir- 
respective of the values of parameters ^k the rational equations (|4.1| )-(|4? 
imply: 



N ~. 

j—\ •'^j •'^k 
N 



N 



j=i 



Xk + 1 



Co- 



Co 



N 

j=l Xk Xj -\- 7 j—^ Xk Xj 



(5.1) 
(5.2) 



and the hyperbolic equations ( |4.1(j| )-( |4.18| ) imply: 

N N 

h E CjCoth.{xj — Xk) ~ CjCoth.{xj — + 7) = 1, 



1. 



(5.3) 
(5.4) 



N N 

— c,coth(xfc — 5j + 7) + X! CjCoth(xfc 
i=i i=i 

These equations are very similar to those found in [0], and can be anal- 
ysed along the same lines. Namely, each pair of equations ( ^.1|) , ( ^.2|) or 
( |5.3| ), ( p.4| ) may be considered as a linear system for 2N variables c^, c^ with 
a Cauchy matrix. This system may be explicitly solved, which allows then 
to write down " Newtonian" equations of motion in terms of x variables only. 
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To appreciate the difference between the equations found in |T3[ and 
here, it will be instructive to compare the Lax representations for them. Our 
results immediately imply the Lax representation of the form 

LM =ML 

where in the rational case 

^7 ^ h 

-f' = X! — ~ T — ^kEjk-, = XI '^^^ — CkEjk-, 

j,k=l ~ +1 j^k=l ~ "^^ 

and in the hyperbolic case 



sinh(7)e^* ^ . . he^'' ^ 

L= • 1 . -^CkEjk, M= —TT^ TCkEjk. 

smh(xj - Xfc + 7) ^.^^ smh(xj - Xk) 

So, they belong to the same hierarchies as their continuous time counterparts; 
the interpolating flows are given by the same formulas as in Theorem 2. 

The Lax representations for the systems found in have the same (up 
to a simple gauge transformation) Lax matrices L, but different matrices M.^ 
namely in the rational case 

7 

M= ^z— —CkEjk, 

j,k=i "^i ^fc + 7 

and in the hyperbolic case 

smh(xj - Xfc + 7) 

Hence, they also belong to the same hierarchies as their continuous time 
counterparts, but are interpolated by different flows. In fact, interpolating 
flows were found in |jl3| to correspond to 

f(L) = -iL + h^~')-' 
in the rational case, and to 

/ L + /i(coth(7)-l)/ \ 
' 2h °^l^L + /i(coth(7) + l)Jy' 
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in the hyperbolic case. 

We see that the maps found in the present Letter serve as approximations 
to the flow of the RS hierarchy with f{L) = L, which corresponds to the 
Hamiltonian function <f{L) = tr(L), while the maps found in 



serve as 



approximations to another flow of the RS hierarchy, characterized by f{L) 
— L~^, which corresponds to the Hamiltonian function ip{L) = tr(L~^). [] 



6 Conclusion 

This paper contains the application of a general procedure of integrable dis- 
cretizations to the rational and hyperbolic spin RS models. There arise 
several natural questions. 

• What about elliptic spin RS model? We have not pursued this point, 
basically because of the lacking Hamiltonian structure of this model. 
Indeed, a natural extension of our results to the elliptic case leads to a 
discrete Lax equation with a spectral parameter A: 

L{X)M{X) =M{X)L{X) 

with the same Lax matrix 

N 

= H ^{xj - Xk + 'j,X){bj\ak)Ejk 
j,k=i 

as for the continuous system, and 

N 

M{X) = ^ - Xk, X){bj\ak)Ejk 
j,k=i 

(here A) = a^x + X)/a{x)a{X), and a{x) is the Weierstrass sigma 
function). However, in absence of the underlying Hamiltonian structure 
we can not make any assertions about complete integrability of the 
resulting system. A surrogate for this - an explicit solution of arbitrary 
flow in a hierarchy - is also not available in the elliptic case. We intend 
to return to the elliptic case in the future. 

^Both maps are derived in the recent preprint along with their Lax representa- 
tions (in the general elliptic case). However, Hamiltonian interpretation and relations to 
continuous time hierarchies are not discussed there. 
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• What about spin CM models ? These models, at least their rational and 
hyperbolic versions, can be discretized along the same lines as the spin 
RS models. However, this (paradoxically) turns out to be a technically 
more difficult problem. A reference to the standard implicit function 
theorem when discussing the equations of the type ( [4.1|) -( ^l3D has to be 
replaced by less trivial algebraic considerations. This will be discussed 
elsewhere. 
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